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Abstract 

Background: Hypermethylation of the suppressor of cytokine signaling 3(SOCS3) promoter has been reported to predict a 
poor prognosis In several cancers including glloblstoma multiforme (GBM). We explored the function of S0CS3 promoter 
hypermethylation In GBM cohorts, including analysis of the CpG Island methylator phenotype (CIMP), when a large number 
of gene loci are simultaneously hypermethylated. 

Methods: k whole genome promoter methylatlon profile was performed In a cohort of 33 GBM samples, with 13 long-term 
survivors (LTS; overall survival > 18 months) and 20 short-term survivors (STS; overall survival < 9 months). The S0CS3 
promoter methylation status was compared between the two groups. In addition, we Investigated the relationship of 
S0CS3 promoter methylation and G-CIMP status. 

/feii//f5/ Interestingly, in our present study, we found that S0CS3 promoter methylation was statistically significantly higher 
In the 13 LTS than that in the 20 STS. Furthermore, high S0CS3 promoter methylation detected via pyro-sequencing 
predicted a better prognosis in an independent cohort containing 62 GBM patients. This correlation was validated by the 
dataset from the Cancer Genome Atlas{TCGA) and the Chinese Cancer Genome Atlas{CGGA). In addition, we found that 
hypermethylation of the S0CS3 promoter was tightly associated with the G-CIMP-posltlve GBM patients. 

Conclusions: Using a total of 359 clinical samples, we demonstrate that S0CS3 promoter hypermethylation status has a 
favorable prognostic value in GBM patients because of whole genome methylation status. Particularly, the 
hypermethylation of the S0CS3 promoter Indicates positive G-CIMP status. 
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introduction 

Glioblastoma is the most malignant primary brain tumor in 
adults with an overall survival rate of about 1.5 years even when 
treated with radical regimens including surgical resection, and 
radiotherapy with concomitant and/or adjuvant temozolomide 
chemotherapy [1]. Although the exact mechanism of GBM 
development and progression is still unknown, certain molecular 
biomarkers are related to tumorigenesis and progression of GBM 
at the genetic, epigenetic, and transcriptional levels [2,3,4,5]. 

However, markers for GBM that have prognostic value in 
signaling transduction pathways have not been fuUy elucidated 
yet.The Janus kinase/signal transducer and activator of transcrip- 
tion 3 (JAK/STAT3) signaling pathways transmits extracellular 
signals into the nucleus where it regulates DNA transcription and 
activity in the cell [6]. The suppressor of cytokine signaling 3 



(SOCS3) is an endogenous inhibitor of the JAK/STAT3 signaling 
pathway, modulating cell activities via suppressing transcription. 
Recently, some studies have reported that SOCS3 functions as a 
tumor suppressor in multiple tumor types, including GBM 
[7,8,9,10]. 

DNA methylation is a precisely regulated process in normal cells 
that becomes drastically modified in cancer cells [5, 11, 12]. 
Hypomethylation of oncogene promoters and hypermethylation 
of tumor suppressor gene promoters are pivotal alterations in 
cancer development [13, 14]. Moreover, DNA methylation is 
typically a stable and inheritable epigenetic pattern that can 
persist for several cell generations, which potentially broadens its 
clinical practical applicability [15]. 

Hypermethylation of oncogenic genes is a favorable indictor for 
GBM patients. A variety of studies have reported that hypermetliylation 
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of the SOCS3 promoter predicts poor prognosis in certain 
cancers, including GBM[16,17,18,19]. However, in our study, 
hypermethylation of the SOCS3 promoter was associated with 
better outcomes for GBM patients. In addition, we found that 
hypermethylation of the SOCS3 promoter in GBM was tightly 
associated with the G-CIMP-positive GBM patients. 

Materials and Methods 

Patients and samples 

All patients with primary GBM were from the Chinese Glioma 
Genome Atlas (CGGA) who underwent surgical resection between 
January 2006 and December 2010 and subsequently received 
radiotherapy and/or adjuvant temozolomide. Tumor tissue 
samples were obtained by surgical resection before the treatment 
with radiation and/ or chemotherapy. Specimens were snap-frozen 
in liquid nitrogen until nucleic acid extraction. We invited two 
independent neuropathologists to evaluate the specimens histo- 
logically. Primary and secondary glioblastoma were distinguished 
based on patients' clinical history record. Written informed 
consents were obtained from the patients (or their families). No 
minors/children patients were included in our research. This 
study was approved by the Ethics Committee of Capital Medical 
University, Beijing, China. 

DNA extraction 

A hematoxylin and eosin-stained frozen section was prepared 
for assessment of the percentage of tumor cells before DNA 
extraction. Only samples with greater than 80% tumor cells were 
selected. Genomic DNA was isolated from frozen tumor tissues 
using the QIAamp DNA Mini Kit (Qiagen) according to the 
manufacturer's protocol. DNA concentration and quality were 
measured using the NanoDrop ND-1000 spectrophotometer 
(NanoDrop Technologies, Houston, TX). 

Genome-wide DNA metliylation profiling 

We used the lUumina Infinium HumanMethylation27 Bead- 
Chip (lUumina Inc.) [20]. The BeadChip contains 27,578 highly 
informative CpG sites covering more than 14,000 human RefSeq 
genes, and allows researchers to investigate all of these sites per 
sample at a single nucleotide resolution. Bisulfite' modification of 
DNA, chip processing and data analysis were performed following 
the manufacturer's manual at the Wellcome Trust Centre for 
Human Genetics Genomics Lab in Oxford, UK. The array results 
were analyzed with the BeadStudio software (lUumina). We have 
deposited our dataset on Gene Expression Omnibus (GEO) and 
the GEO accession number is GSE53228. 

Pyrosequencing analysis of S0CS3 

Pyrosequencing was supported by Genetech (Shanghai, China) 
and performed using the PyroMark Q96 ID System (Qiagen) 
according to the manufacturer's protocol. Bisulfite modification of 
the DNA was accomplished using the EpiTect Kit (Qiagen). The 
beta value is a quantitative measure of DNA methylation levels of 
specific CpGs using the ratio of intensities between methylated and 
unmethylated alleles[21]. 

Statistical analysis 

T-tests were performed using GraphPad Prism 5. Kaplan- 
Meier survival curves were obtained, and differences in the overall 
survival were tested for statistical significance using the log-rank 
test (GraphPad Prism 5). P<0.05 was considered significant. 



Results 

Hypermethylation of the S0CS3 promoter predicts better 

prognosis for GBM patients 

Our test cohort consisted of 1 3 long-term survivors (LTS) whose 
overall survivals were more than 18 months, and 20 short-term 
survivors (STS) whose overall survivals were less than 9 months. 
The (3 values of the SOCS3 promoter of the two groups are shown 
in Figure lA, which displayed statistically significant difference 
(P<0.01) between the STS group and the LTS group. Thus we 
inferred that hypermethylation of the SOCS3 promoter may 
correlate with favorable prognosis in GBM patients. Our results 
were validated in an independent cohort containing 62 GBM 
samples from Tiantan Hospital, Beijing, China. According to 
average methylation values measured by pyrosequencing (Figure 
SI), 62 samples of the independent validation cohort were divided 
into three groups (Figure IB): average methylation values <30%, 
average methylation values between 30% to 60%, and average 
methylation values > 60%. A comparison of the three groups 
demonstrated that there was a statistically significant difference in 
the patients' survival (P = 0.04) among the groups. Our findings 
were validated using the Cancer Genome Atlas (TCGA) dataset 
(n = 264). In the TCGA validating cohort, we divided all the 
samples into five groups according to their P values. A Kaplan- 
Meier curve of the survival of these 264 patients is shown in 
Figure IC. Significant difference was found between the groups 
with a P value>80% and the groups with relatively lower P values. 

Hypermethylation of the S0CS3 promoter is associated 
with G-CIMP-positive GBM patients 

In addition, we found that hypermethylation of the SOCS3 
promoter is tightly associated with G-CIMP-positive GBM 
patients in two independent cohorts, the CGGA GBM and the 
TCGA GBM cohorts. In the CGGA cohort, the P value of the G- 
CIMP-positive group was 0.66, which is significantiy higher than 
0.26, the P value of the G-CIMP-negative group (Figure 2A, P< 
0.01). Similarly, statistically significant difference was observed in 
the G-CIMP-positive group compared with the G-CIMP-negative 
group in the TCGA GBM cohort with p values of 0.81 and 0.41, 
respectively (Figure 2B, P<0.01). In another TCGA cohort of 
which all the 242 GBM samples were G-CIMP-negative, we 
divided the samples into five groups according the average P 
values, similar to the TCGA validating cohort (Figure IC). Our 
analysis demonstrated that the P value provided little clinical 
prognostic value for these patients as shown in the Kaplan-Meier 
curve (Figure 2C, P = 0.6). These results indicate that the 
prognostic value of hypermethylation of the SOCS3 promoter 
was tightiy associated only with G-CIMP-positive GBM samples. 

A novel indicator for G-CIMP-positive GBM patients 

We performed a receiver operating characteristic curve (ROC 
curve) between the hypermethylation of the SOCS3 promoter and 
G-CIMP to define the exact relationship. According to our data 
analysis, statistical significance was observed in CGGA samples 
(AUG = 0.951, P = 0.001) (Figure 3A). These results indicate a 
robust relationship between the h)rpermethylation of the SOCS3 
promoter and G-CIMP positive. Thus, hypermethylation of the 
SOCS3 promoter is a de novo indicator for G-CIMP. To validate 
our results, TCGA samples were subsequently used to verity this 
relationship; the results are even better than that of CGGA 
samples (AUG = 0.943, P<0.001) (Figure 3B). 
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Figure 1. SOCS3 methylation status in STS and LTS group and vaiidiation cohort. A. In the CGGA G6M cohort, the |3 values of the S0CS3 
promoter of the two groups (STS group and LTS group) are significantly different (P<0.01). B. In an independent validation cohort, survival analysis 
showed that three groups divided by average methylation values are significantly different (P<0.04). C. In the TCGA GBM cohort, the group with |3 
value>80 percent (red) has a significantly longer survival than the other four groups (P = 0.02). 
doi:1 0.1 371 /journal.pone.0091 829.g001 



Discussion 

Glioblastoma is the most malignant primary brain tumor in 
adults, with insidious development, rapid progression and poor 
outcomes. Alterations in cell signaling pathways may be associated 
with the development and progression of GBM. Some prognostic 
bio-markers involved in signaling pathways have been identified. 
Hypermethylation of the SOCS3 promoter has been associated 
with a poor outcome for GBM. From our present research, we 
draw an opposite conclusion to previous studies and show that 
hypermethylation of the SOCS3 promoter predicts an improved 
prognosis for GBM patients. 

DNA metiiylation is a common regulatory process which 
influences cell activities including transcription in normal cells. 
DNA methylation frequently becomes drastically aberrantly 
altered in cancer cells [12]. Hypomethylation of oncogene 
promoters and hypermethylation of tumor suppressor gene 
promoters are pivotal alterations in cancer development. 

The CpG island methylator phenotype (GIMP) is a methylation 
status when a large number of gene loci are simultaneously 
hypermethylated, probably as consequence of mutations of 
methyltransferases or histone-modifying proteins[22], aging[22], 
virus exposure[23,24], chronic inflammation [2 5, 2 6] or other 
underlying factors. Reportedly, GIMP was observed in many 
tumors, including colorectal cancer [2 7, 2 8], adrenocortical 



carcinomas [29], gastric tumors[30,31], liver cancer[23], esopha- 
gus cancer[32], ovarian cancers[33] and acute myelogenous 
leukemia[34,35]. In different tumors, GIMP of the whole tumor 
genome affects different specific genes and functions differently, 
either as favorable or unfavorable predictors for patients. Poorer 
outcome was observed in patients who suffered adrenocortical 
carcinomas with the existence of CIMP[29]. Nevertheless, 
according to previous research, in gastric carcinoma, the prognosis 
of the patients without GIMP was significantly worse compared 
with that of patients with GIMP[29]. Such evidence confirms the 
fact that hypermethylation of the whole cancer genome does not 
necessarily mean better or worse outcomes for patients. Instead, it 
is the specific genes that are aberrantly methylated that determine 
outcomes [2 7]. 

G-CIMP is enriched in a subgroup of glioma, the proneural 
subgroup, according to the TGGA classification scheme for 
glioma[36]. In G-GIMP-positive samples, a large number of GpG 
island loci located in specific gene promoters are hypermethylated 
and patients usually have better outcomes [2 7]. According to our 
research and data analysis, hypermethylation of the SOCS3 
promoter is highly associated with G-CIMP-positive samples and 
predicts improved outcomes for patients, but is not a predictor for 
G-GIMP-negative patients. Therefore, we conclude that SOGS3 
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Figure 2. G-CIWIP status in CGGA and TCGA cohort. A. In the CGGA GBM cohort, the (3 values of the S0CS3 promoter of the two groups (G- 
CIMP-positive group and G-CIMP-negative group) displayed statistically significant difference (P<0.01). B. In the TGGA GBM cohort, the (3 values of 
the S0CS3 promoter of the two groups (G-CIMP-positive group and G-CIMP-negative group) also displayed statistically significant difference 
(P<0.01). C. In the G-CIMP-negative TCGA samples, there was no significant difference among the five groups (P = 0.60). 
doi:1 0.1 371/journal.pone.0091 829.g002 
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Figure 3. Hypermethylation of SOCS3 promoter is consistent witii G-CIIMP. A. In the CGGA samples, a ROC curve demonstrated a tight 
relationship between hypermethylation of the S0CS3 promoter and G-CIMP-positive status (AUC = 0.951, P = 0.001). B. In theTCGA samples, a similar 
relationship between hypermethylation of the S0CS3 promoter and G-CIMP-positive status is shown in the ROC curve (AUC = 0.943, P<0.001). 
doi:1 0.1 371/journal.pone.0091 829.g003 



hypermethylation status has favorable prognostic value in GBM 
patients because of its whole genome methylation status. 

SOCS3 functions as a tumor suppressor in many cancers 
including GBM. According to the bio-effects of the genetic 
hypermethylation process, hypermethylation of tumor suppressor 
gene promoters theoretically is aversive for tumorigenesis or 
progression. Furthermore, many studies have confirmed the effect 
of SOCS3 in GBM samples. In G-CIMP-positive samples, as our 
data showed above, the SOCS3 promoter is hypermethylated 
along with a variety of other loci. The hypermethylation of the 
SOCS3 promoter is just a part of the whole genome methylation 
status and its negative effect on tumorigenesis or progression may 
be neutralized by the comprehensive genome hypermethyla- 
tion [3 7]. This hypothesis may explain why hypermethylation of 
the SOCS3 promoter predicts favorable prognosis in GBM 
patients. In addition, other potential signaling pathways may be 
uncovered for which hypermethylation of the SOCS3 promoter 
serves as a better prognosticator. Because this single gene 
alteration accompanies whole genome hypermethylation, SOCS3 
can be regarded as a pivotal gene that functions as a predictor for 
the whole genome methylation status (G-CIMP). As we revealed in 
this research, SOCS3 hypermethylation is a de novo indicator for 
G-CIMP and predicts better patients' outcomes. The prognostic 
value of SOCS3 hypermethylation is also practical as it is easy to 



perform in clinical practice and could be helpful in determining 
therapeutic regimens for GBM patients. 

Conclusions 

In summary, we found that hypermethylation of the SOCS3 
promoter predicts favorable prognosis. Our results were vahdated 
in an independent cohort containing 62 GBM samples as well as in 
a TCGA GBM cohort. Further investigation is needed to uncover 
the exact mechanism of how hypermethylation of SOCS3 
promoter affects the normal processes in the cell and its 
relationship to tumorigenesis and progression. We also found that 
SOCS3 is a de novo indicator for G-CIMP-positive GBM patients. 

Supporting information 

Figure SI Pyrosequencing for SOCS3 promoter meth- 
ylation. This figure shows unmethylated and methylated SOCS3 
promoters using pyrosequencing. 
(TIF) 
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